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ABSTRACT 

Measurements  on  long-period  P,  S,  and  LR  phases  were  extracted  from  an 
eight-month  data  base  for  nine  stations  of  the  Long  Period  Experimental  network. 
Amplitude-distance  curves  were  plotted  for  P and  S from  shallow-focus  events; 
these  curves  generally  agreed  with  the  commonly-used  Gutenberg-Richter  B factors 
The  amplitude-distance  curve  for  LR  wai  different  in  slope  than  the  commonly- 
used  correction  curve  for  Ms;  however,  this  difference  would  not  be  sufficient 
to  affect  routine  network  Ms  estimates  by  more  than  roughly  .1  magnitude  unit. 
The  Q for  20-second  LR  implied  by  this  new  data  is  roughly  700,  in  contrast  to 
the  300  implied  by  the  accepted  amplitude-distance  curve. 
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INTRODUCTION 


A large  data  base  of  visual  long-period  measurements  on  recordings  of  the 
Long  Period  Experimental  (LPE)  network  of  seismometers  has  been  reported  by 
von  Seggern  (1976).  This  network  has  the  advantage  of  high  magnifications  due 
to  rigid  environmental  control  and  emplacement  at  depth  for  the  seismometers 
(Pomeroy  et  al. , 1969).  Thresholds  for  long-period  phases  were  sufficiently 
low  that  copious  data  was  accumulated  at  the  SDAC  while  analyzing  the  first 
eight  months  of  recordings  in  1972  for  the  nine  LPE  sites  then  operating; 
these  LPE  thresholds  were  (in  terms  of  90%  detection  by  >_  2 stations)  roughly 
5.7  m^  for  P,  5.5  for  S,  5.4  for  LQ,  and  5.0  for  LR,  considering 
detection  of  all  events  worldwide  (von  Seggern,  1976).  The  nine  stations 
operating  during  this  period  were:  ALQ  (Albuquerque,  New  Mexico),  CHG  (Chiang 

Mai,  Thailand),  CTA  (Charters  Towers,  Australia),  EIL  (Eilot,  Israel),  FBK 
(Fairbanks,  Alaska),  KIP  (Kipapa,  Hawaii),  KON  (Kongsberg,  Norway),  OGD 
(Ogdensberg,  New  Jersey),  and  TLO  (Toledo,  Spain).  Station  KIP  began  operation 
on  March  6,  1972;  and  FBK  operated  only  until  April  26,  1972. 

Drawing  on  that  LPE  data  base,  this  report  will  present  the  amplitude- 
distance  relation  of  the  two  major  long-period  body-wave  phases,  P and  S,  and 
of  LR  near  a period  of  twenty  seconds.  The  LQ  data  was  not  considered  sufficient 
in  the  period  range  near  20  seconds  where  Mg  is  usually  obtained.  This  deficiency 
occurred  because  available  horizontal  recordings  were  seldom  aligned  along  the 
radial  and  transverse  directions  of  motion,  thus  contaminating  shorter  periods 
of  LQ  motion  with  Rayleigh  waves.  In  addition  noise  on  the  horizontal  seismo- 
graph recordings  was  usually  higher  than  on  vertical  ones.  The  data  base  of 
long-period  measurements  was  associated  with  a merged  list  of  NEIS,  LASA 
bulletir , and  NORSAR  bulletin  epicenters.  Collation  of  presumably  identical 


von  Seggern,  D.  H. , 1976,  Final  report  on  the  analysis  of  recordings  from  the 
Very  Long  Period  Experimental  stations,  SDAC-TR-76-1,  Teledyi.e  Geotech, 
Alexandria,  Virginia  (in  press). 

Pomeroy,  P.  W.,  G.  Hade,  J.  Savino,  and  R.  Chander,  1969,  Preliminary  results 
from  high-gain  wide-band  long-period  electromagentic  seismograph  systems, 
J.  Geophys.  Res.  v.  74,  p.  3295-3298. 


epicenters  was  accomplished  using  criteria  of  less  than  10°  spatial  separation 
and  less  than  one  minute  time  separation.  No  depth  criterion  was  used  in 
collating  because  of  the  general  lack  of  such  information  in  the  LASA  and 
NORSAR  bulletins.  The  merged  list  comprised  5390  epicenters  for  the  first 
eight  months  of  1972.  Over  one-half  of  these  epicenters  had  phases  recorded  at 
one  or  more  of  the  LPE  stations. 


BODY  WAVES 


Procedure 

In  the  LPE  data  base  (von  Seggern,  1974)  several  thousand  P amplitudes 
and  S amplitudes  were  available  for  inferring  long-period  amplitude-distance 
relations.  Only  those  phases  from  shallow-focus  (<  60  km  earthquakes  are 
considered  because  the  data  from  deeper  events  was  not  considered  sufficient 
to  establish  depth  corrections.  Since  these  phases  represented  events  whose 
size  extended  over  two  or  three  orders  of  magnitude,  some  manner  of  normalizing 
the  body-wave  amplitudes  among  events  was  necessary  to  obtain  meaningful 
displays  of  amplitude  versus  distance  for  these  phases.  We  chose  to  normalize 
the  body-wave  amplitudes  by  the  average  Mg  for  the  event,  provided  at  least 
three  stations  had  an  Mg  estimate  for  the  event.  Body  waves  for  events  which 
did  not  meet  this  criteria  were  discarded.  However,  only  a small  part  of  the 
available  P and  S amplitudes  were  eliminated  since  the  long-period,  body-wave 
thresholds  were  at  least  one-half  magnitude  unit  higher  than  the  LR  threshold 
(von  Seggern,  1976),  and  since  almost  all  events  which  had  one  or  more  reported 
long-period  body  waves  also  had  at  least  three  stations  detecting  LR.  the 
level  of  LR  generation  should  correlate  well  with  that  of  the  long-period  body- 
wave  generation  by  a shallow  earthquake  according  to  seismic  source  theory 
(Blandford,  1975;  Gilbert,  1973;  Douglas  et  al. , 1973).  Thus,  normalizing 
long-period  body-wave  amplitude  by  M (log  A/T-M  ) should  remove  most  of  the 
scatter  due  to  source  factors.  The  number  of  amplitudes  able  to  be  normalized 
in  this  manner  were  1279  for  P waves  and  2226  for  S wa.'es.  No  restriction  was 
placed  on  period;  the  range  extended  from  5 to  80  seconds,  with  the  majority  of 
measured  periods  falling  in  the  20-30  second  range,  which  is  considerably  longer 
than  typical  periods  for  long-period  body  waves  seen  on  WWSSN  recordings. 


Blandford,  R.  R. , 1975,  A source  theory  for  complex  earthquakes.  Bull.  Seism. 
Soc.  Amer. , in  press. 

Gilbert,  F.,  1973,  The  relative  efficiency  of  earthquakes  and  explosioi  s in 
exciting  surface  waves  and  body  waves,  Geophys.  J.  R.  Astr.  Soc..  v.  33 
p.  487-488. 

Douglas,  A,,  J.  A.  Hudson,  and  C.  Blarney,  1973,  A quantitative  evaluation 

of  seismic  signals  at  teleseismic  distance— III  computed  P and  Rayleigh 
wave  seismograms,  Geophys.  J.  R.  Astr.  Soc.,  v.  28,  p.  385-410. 
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Results 


A plot  of  these  amplitudes  versus  distance  is  shown  in  Figure  1 for  P 
waves  and  Figure  2 for  S waves.  Amplitudes  were  grouped  into  2°  bins  and 
averaged;  the  standard  deviation  of  the  mean  amplitudes  are  indicated  by  error 
bars.  For  comparison  we  have  included  the  Gutenberg-Richter  (1956)  tabulated 
amplitude  terms  for  P and  S waves.  We  note  that  their  curves  were  compiled 
with  phases  having  shorter  periods  than  those  in  our  data  base;  although  the 
period  range  is  not  specifically  stated,  Gutenberg  (1945a)  implies  that  data  for 
these  curves  had  periods  generally  between  1 and  10  seconds. 

Considering  now  the  P-wave  curves  in  Figure  1,  we  note  that,  although  the 
general  trend  is  similar,  significant  departures  in  our  data  from  Gutenberg  and 
Richter's  occurs  in  the  ranges  15°-18°,  38°-44°,  and  91°-99°.  For  15°-18°  an 
explanation  can  probably  be  found  in  regional  effects  on  amplitudes  as  the 
extent  of  the  low  velocity  zone  is  critical  here.  No  data  is  shown  for  A < 15° 
because  of  large  standard  deviations.  For  38°-44°,  there  is  no  satisfactory 
explanation  of  the  high  values  obtained  by  Gutenberg  and  Richter  since  geo- 
metrical spreading,  the  most  important  influence  for  long  period  P waves, 
causes  only  a very  smooth,  almost  neglibible,  decay  in  this  range  (Booth  et  al., 
1974).  For  91  -99°,  our  decay  is  not  as  great  as  that  suggested  by  Gutenberg 
end  Richter;  recalling  that  our  periods  were  quite  long  compared  to  Gutenberg 
a.A  Richter,  we  can  explain  this  facet  of  our  data  by  the  diffraction  theory 
of  Phinney  and  Cathles  (1969)  for  ? waves  near  the  core.  Overall,  our  curve 
is  more  similar  to  those  of  Willey  et  al.  (1970)  and  Booth  et  al.  (1974). 


Gutenberg,  B. , and  C.  F.  Richter,  1956,  Magnitude  and  energy  of  earthquakes, 
Annali  de  Geofisica  r 9,  p.  1-15. 

Gutenberg,  B. , 1945a,  Amplitudes  of  P,  PP,  and  S and  magnitude  of  shallow 
earthquakes.  Bull.  Seism.  Soc.  Am.,  v.  35,  p.  57-69. 

Booth,  D.  C.,  P.  D.  Marshall,  and  J.  B.  Young,  1974,  Long  and  short  period 
P-wave  amplitudes  from  earthquakes  in  the  range  0°-114°,  Geophys.  J.  R. 
Astr.  Soc.,  v.  39,  p.  523-537. 

Phinney,  R.  A.,  and  L.  M.  Cathles,  1969,  Diffraction  of  P by  the  core:  a study 
of  long-period  amplitudes  near  the  edge  of  the  shadow,  J.  Geophys.  Res., 
v.  74,  p.  1556-1574. 

Willey,  G. , J.  R.  Cleary,  and  P.  D.  Marshall,  1970,  Comparison  of  least  squares 
analysis  of  long  and  short  period  P wave  amplitude,  Geophys.  J.  R.  Astr. 
Soc.,  v.  19,  p.  439-445. 
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Turning  now  to  the  S-wave  curves  in  Figure  2,  we  see  that  our  data  agrees 
with  the  Gutenberg-Richter  curve  satisfactorily  except  for  A < 25°  where,  as 
tor  P waves,  the  regional  affects  of  the  low-velocity  zone  probably  are  causing 
the  divergence.  Data  for  A < 13  is  not  plotted  because  of  large  standard 
deviations.  Again,  there  are  no  theoretical  reasons  for  fluctuations  in  the 
curves  beyond  25  and  in  fact,  for  our  data,  a very  smooth  curve  could  be  drawn 
to  lie  within  the  95%  confidence  limits  of  all  the  points  (tnese  limits  would 
be  approximately  twice  the  standard  deviation  intervals  shown).  Other  than  the 
Gucenberg-Richter  curve  shown,  there  is  no  S-wave  amplitude  curve  over  this 
distance  range  to  use  as  a basis  of  comparison. 


RAYLEIGH  WAVES 


Data 

For  this  section  on  LR,  we  shall  describe  the  data  base  in  more  detail. 

The  LR  measurements  used  in  this  study  were  made  peak- to- trough  on  the  maximum 
of  the  Rayleigh-wave  group  in  a period  of  17-23  seconds.  Roughly  two-thirds 
of  the  data  fell  in  the  range  19-21  seconds.  Many  measurements  were  also  made 
on  LR  outside  these  limits  but  are  ignored  in  this  study.  The  visually 
measured  amplitudes  were  reduced  to  ground  displacement  using  available  magni- 
fication curves,  and  unlike  the  P and  S amplitudes,  were  not  divided  by  period 
for  our  purposes.  We  accepted  measurements  on  horizontal  recordings  if  none 
were  available  from  the  vertical,  but  only  a tiny  fraction  of  the  LR  data  fell 
in  this  category.  We  also  restricted  our  data  to  that  associated  with  shallow- 
focus  (<_  60  km)  earthquakes;  this,  however,  eliminated  only  a small  part  of  the 
available  data. 

The  data  in  general  was  considered  to  be  of  high  quality  since  it  was 
assembled  by  experienced  analysts  who  demonstrated  a very  low  false  alarm  rate 
(von  Seggern,  1974)  and  who  were  operating  in  response  to  predicted  LR  arrival 
times  from  known  epicenters  made  with  a global  group-velocity  grid  (Kimball, 
1969).  It  might  be  expected  though  that  amplitudes  near  the  threshold  of 
detection  were  contaminated  by  noise  or  were  possible  false  alarms;  such  ampli- 
tudes would  be  reported  for  the  smaller  events  and  for  larger  epic">ntral 
distances.  We  were  in  fact  forced  by  our  choice  of  statistical  analysis,  as 
discussed  in  the  next  section,  to  limit  the  data  to  those  events  which  had  at 
least  five  LR  amplitudes,  restricted  as  stated  above,  in  order  to  remain  within 
the  computer  core  capacity.  This  alone  eliminated  roughly  one-half  of  the 
available  data;  but  those  events  with  four  or  less  detections,  being  of  smaller 
magnitude,  would  include  most  of  the  reported  amplitudes  of  less  quality. 

There  remained  4456  LR  amplitudes  from  691  events  for  the  analysis  after  all 
the  criteria  which  have  been  mentioned  were  applied  in  restricting  the  data 
base. 


Kimball,  B.  C.,  1969,  Prediction  of  seismic  surface-wave  travel  times,  Technical 
Report  69-40,  Teledyne  Geotech,  Garland,  Texas. 
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Statistical  Model 


We  have  followed  Carpenter  et  al.  (1967)  in  using  the  model: 

Am  ' si  + Fj  + \ + eijk  a) 

where 

= logarithm  of  observed  amplitude, 

S±  = station  effect, 

Fj  = source  effect, 

R^,  = distance  effect, 

6ijk  = random>  normally  distributed,  error. 

One  >can  solve  for  the  effeccs  and  their  variances  by  least  squares  if  the  side 
conditions  = 0 and  = 0 are  imposed  (Scheffe,  1959). 

Hie  results  from  applying  the  above  model  must  be  judged  against  the 
adequacy  of  the  model  itself.  Note  that  we  have  ignored  any  interactions  which 
might  arise  due  to  non-unifcrmity  of  sources  (radiation  patterns  from  double 
couples)  or  from  the  heterogeneity  of  the  earth  (path  effects  on  amplitude). 
Ignoring  such  interactions  would  be  especially  worrisome  if  the  data  were  not 
evenly  distributed  as  regards  source,  station,  and  distance.  Tables  I and  II 
show  the  distribution  of  the  data  used  in  the  least-squares  solution  for 

the  effects  in  equation  (1) ; examination  of  these  tables  reveals  that  no  great 
concentrations  of  data  occur,  and  therefore  it  seems  apparent  no  large  biases 
in  the  results  should  occur.  We  do  point  out  some  irregularities  though,  such 
as  the  large  proportion  of  ALQ,  OGD,  and  TLO  data  at  8G°  to  120°  and  a similar 
happenstance  for  CHG  and  CTA  at  10°  to  60°.  These  observed  concentrations  are 
of  course  due  to  the  non-random  pattern  of  seismicity  over  the  earth.  If 
significant  interactions  of  opposite  sign  occurred  between  station  and  distance 
for  these  two  groups,  then  the  overall  trend  of  the  distance  effects  could 
be  seriously  biased.  We  do  not,  however,  have  sufficient  geophysical  knowledge 


Carpenter,  E.  W. , P.  D.  Marshall,  and  A.  Douglas,  1967,  The  amplitude  distance 
curve  for  short-period  teleseismic  P waves,  Geophys.  J.  R.  Astr.  Soc 
v.  13,  p.  61-70.  * 

Scheffe,  H. , 1959,  The  Analysis  of  Variance.  John  Wiley  and  Sons,  Inc.,  New 
York,  New  York. 


-14- 


TABLE  I 


Distribution  or  Observations  by  Distance  and  Station 


Distance 

ALQ 

CHG 

CTA 

Number  of 
EIL  FBK 

0°-10° 

3 

2 

0 

3 

2 

10°-20° 

10 

7 

63 

24 

8 

20°-30° 

25 

61 

40 

9 

17 

30 °-4Q° 

28 

35 

72 

17 

10 

40 °-50° 

20 

66 

48 

15 

17 

50 °-60° 

70 

48 

29 

32 

18 

60°-70° 

32 

20 

25 

32 

29 

O 

O 

00 

1 

o 

o 

37 

35 

40 

46 

39 

o 

O 

ON 

1 

O 

o 

00 

69 

18 

39 

85 
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to  auantify  the  magnitude  of  such  interactions,  and  the  inclusion  of  them  in 
the  model  would  make  it  intractable  for  solution  with  our  large  body  of  data. 

Results 

The  station  effects  are  estimated  using  equation  (1)  and  are  listed  in 
Table  III;  the  confidence  limits,  not  listed,  were  all  on  the  order  of  +.02 
from  the  mean  at  the  95%  level.  The  negative  of  these  quantities  would  be 
applied  as  a correction  to  calculation  of  Mg  for  routine  seismol ogical  purposes 
although,  except  for  OGD,  the  values  are  so  small  that  applications  of  these 
corrections  to  network  estimates  of  Mg  would  seem  superfluous,  especially 
since  other  effects,  such  as  due  to  path  and  source  non-uniformities  and  to 
measurement  error,  will  dominate.  The  station  effects  of  Table  III  reflect  not 
just  the  structure  near  each  station  but  the  overall  path  effects  for  the 
signals  recorded  at  each  one,  and  so  we  cannot  readily  correlate  the  observed 
effects  with  those  effects  on  Rayleigh-wave  amplitudes  predicted  from  the 
crustal  structure  at  each  site. 

The  solution  for  the  R^  effects  in  equation  (1)  are  illustrated  in 
Figure  3.  The  average  global  diminution  of  20-second  LR  with  distance  is 
represented  by  these  terms.  We  show,  for  comparison,  the  -1.66*logA  relation, 
arbitrarily  scaled  along  the  ordinate,  which  is  the  accepted  (Vanek,  1962) 
relation  for  normalizing  LR  amplitudes  in  routine  calculations  of  Mg  according 
to: 

= log  A + 1.66’logA  - 130 

where  A is  in  millimicrons  peak-to-trough  and  A is  in  degrees.  The  amplitude 
relation  was  first  proposed  by  Gutenberg  (1945b)  in  a slightly  different  form 
for  15°  < A < 140°.  Quite  clearly,  the  LPE  data  deviates  significantly  from 
this  line;  and  we  have  fitted  the  unweighted  terms  from  15°  to  135°  inclusive 
with  a least-squares  line  as  shown.  The  95%  confidence  limits  on  the  slope 
were  1.08  + .04.  A proposed  Mg  formula,  based  on  the  LPE  data,  is  given  by 

Mg  = log  A + 1.08* logA  - 0.22 

Vanek,  J.,  A.  Zatopek,  V.  Karnik,  N.  V.  Kondorskaya,  Yu.  V.  Riznichenko,  E.  R. 
Savarensky,  S.  L.  Solov'ev,  and  N.  V.  Shebalin,  1962,  Bulletin  of  the 
Academy  of  Sciences,  USSR,  Geophysics  Series,  No.  2 (February  1962), 
p.  108-111  (English  translation). 

Gutenberg,  B.,  1945b,  Amplitudes  of  surface  waves  and  magnitudes  of  shallow 
earthquakes.  Bull.  Seism.  Soc.  Amer.,  v.  35,  p.  3—12. 


The  constant  was  set  by  minimizing  the  average  difference  between  old  and  new 

Mg  estimates  for  the  691  events  in  our  data  set;  he  two  formulas  give  identical 

M at  A = 72°. 
s 

The  R^  terms  from  15°  to  175"  can  be  fit  by  the  theoretical  relation  for 
surface  waves  propagating  on  a sphere;  this  relation  involves  the  stationary- 
phase  approximation,  geometrical  spreading,  and  attenuation  (Sato,  1967): 


-?rr 

R - „v-1/2  • _1/2A  QUT 
K,  - r r sm  Ae 

where  c = constant 

r = epicentral  distance  in  km, 

A = epicentral  distance  in  degrees, 

Q = quality  factor, 

T = period, 

U = group  velocity. 

We  shall  set  T = 20  seconds  and  U = 3.5  km/sec,  a value  appropriate  as  a global 
average  for  propagation  of  LR  waves  at  a period  of  twenty  seconds.  Taking 
logarithms  and  rearranging  the  above,  we  have 

Rk  + 2^  log  r + y log  (sinA)  = log  c + (^)  (-.0195r). 

This  can  be  expressed  in  the  simple  linear  regression  form: 
yk  - a + Hxk 

in  order  to  solve  for  Q = 1/6.  A fit  to  the  data,  as  shown  in  Figure  3, 
assuming  no  error  in  the  R^  terras  and  applying  equal  weights,  resulted  in 
Q * 704  with  95%  confidence  limits  of  597  and  857. 

Discussion 

There  is  a serious  discrepancy  in  the  rate  of  amplitude  diminution  for 
20-second  LR  between  Gutenberg's  (1945b)  data  and  the  LPE  data.  There  is  no 
other  published  study  on  the  scale  of  Gutenberg's  with  which  to  compare  these 
present  results.  We  point  out  though  that  von  Seggern  (1975)  has  determined 


Sato,  R. , 1967,  Attenuation  of  seismic  waves,  J.  Phys.  Earth,  v.  15,  p.  32-61. 

von  St-ggern,  D.  H.,  1975,  Q for  twenty-second  Rayleigh  waves  from  complete 

great-circle  paths,  SDAC-TR-75-3,  Teledyne  Geotech,  Alexandria,  Virginia. 
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a Q for  20-second  LR  on  the  order  of  500  using  many  great  circle  paths  around 

the  globe,  and  such  an  estimate  of  Q would  agree  with  a value  roughly  halfway 

between  those  implied  by  the  LPE  data  and  by  the  Gutenberg  data.  It  is 

impossible  to  determine  if  any  significant  bias  exists  in  the  result  of 

Gutenberg  since  he  does  not  list  the  epicenters  used  o.'  the  distribution  of 

observation  distances  for  each  station;  but  with  a group  of  33  world-wide 

stations  and  over  1000  observations,  it  would  seem  that  the  data  is  sufficient 

to  suppress  any  serious  bias  in  his  result.  By  contrast,  we  used  9 stations  ^ 

and  over  4000  observations,  again  seemingly  sufficient  to  suppress  any  serious 

bias.  In  both  Gutenberg's  and  the  present  study,  only  periods  near  20  seconds 

were  allowed;  and  although  Gutenberg  took  measurements  from  the  horizontal 

recordings,  this  should  ma’.;°  no  difference.  We  postulate  then  that  the  true 

amplitude  diminution  lies  somewhere  in  between  and  that  unidentifiable  station- 

path-source  interactions  peculiar  to  each  study  has  caused  the  results  to 

deviate  in  opposite  directions  from  the  underlying,  true  relation  for  LR 

diminution  with  distance.  We  emphasize  that  the  differences,  in  the  context 

of  M determination,  are  probably  insignificant;  for  as  Figure  3 shows,  the 
s 

maximum  difference  between  the  -1.66  log  A and  -1.08  log  A relation  would  only 
be  about  .1  magnitude  unit  so  that,  in  routine  network  estimation  of  Mg,  values 
calculated  by  both  formulas  would  agree  very  closely  in  all  cases. 
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CONCLUSIONS 


Amplitude-distance  curves  formed  using  long-period  P and  S waves  from  a 
large  LPE  data  base  generally  agreed  with  those  published  by  Gutenberg  and 
Richter  (1956"> . In  the  case  of  long-period  P,  the  LPE  data  agreed  more  with 
recent  studies  using  LRSM  and  WWSSN  data.  In  the  case  of  long-period  S,  no 
other  satisfactory  curve  f jr  comparison  is  available.  But  in  neither  case 
are  the  curves  given  ’ere  or  elsewhere  sufficiently  well-determined  to  promote 
abandoning  the  Gutenberg-Richter  curves  for  a newer  standard.  Some  smoothing 
of  their  curves  based  on  geometrical  spreading  arguments  might  be  desirable 
but  still  would  not  significantly  affect  network  magnitude  determinations 
using  these  phases.  Our  data  appears  to  support  the  core  diffraction  affects 
on  amplitude  proposed  theoretically  by  Phinney  and  Cathles  (1969). 

The  amplitude-distance  curve  for  20-second  LR  derived  from  the  LPE  data 
is  significantly  different  from  Gutenberg's  (1945b),  amounting  to  +.15  M 

s 

unit  at  15°  and  sloping  down  to  -.09  Mg  unit  at  ltD°.  There  is  no  obvious 
reason  for  this  discrepancy.  However,  in  the  context  of  network  estimation 
of  Mg,  this  discrepancy  can  be  safely  ignored,  with  either  Gutenberg's  or  the 
LPE  formula  giving  nearly  equivalent  average  M in  all  cases.  The  two  formulas 
imply  widely  differing  Q on  a global  scale  though,  roughl’  300  for  Gutenberg's 
and  700  for  the  LPE  one;  thus  we  feel  that  the  two  sets  ot  data  represent 
sampling  of  the  earth's  crust  which  cannot  overlap  to  a great  degree.  A con- 
tribution to  resolving  this  discrepancy  could  only  come  from  a study  using  a 
large  body  of  world-wide  data;  ouch  data  could  be  taken  from  the  present  NEIS 
files  or  possibly  from  the  future  files  created  by  the  Network  Event  Processor 
at  the  SDAC. 
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